多 变量 追踪 研究 的 模型 整合 与 拓展 : 
考察 往复 式 影 响 与 增长 趋势 
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Ж 要 追踪 研究 当中 ， 交 又 滞后 模型 可 以 探 完 多 交 量 之 间 往 复式 影响 ， 潜 增长 模型 可 以 探 
究 个 体 增 长 趋势 。 对 两 类 模型 进行 整合 ， 例 如 同时 关注 往复 式 影响 与 个 体 增长 趋势 ， 同 时 可 
以 定义 测量 误差 、 随 机 截 距 等 变异 成 分 ,衍生 出 随机 蕉 距 交 又 灌 后 模型 、 特 质 -状态 -误差 模 
型 、 自 回归 潜 增 长 模型 、 结 构 化 残 差 潜 增 长 模型 等 。 以 交叉 滞后 模型 和 潜 增 长 模型 分 别 作为 
基础 模型 ， 从 个 体 间 /个 体内 变异 分 解 的 角度 对 上 述 各 类 模型 杭 理 ， 整 合 出 此 类 模型 的 分 析 框 
架 ， 并 拓展 建立 “因子 结构 化 潜 增 长 模型 (factor latent curve model with structured reciprocals)” 
作为 统合 框架 。 通 过 实证 研究 (早期 儿童 的 追踪 研究 -幼儿 园 版 ,ECLS-K), 建立 21049 名 儿童 
的 阅读 和 数学 能 力 的 往复 式 影响 与 增长 趋势 。 研 究 发 现 ， 分 离 了 稳定 特质 的 模型 拟 合 最 优 。 
研究 也 对 模型 建 模 思 路 和 模型 选择 提供 了 建议 。 

关键 词 追踪 研究 ， 往 复式 影响 ， 增 长 趋势 ， 因 子 结构 化 潜 增 长 模型 


分 类 号 B841 


追踪 研究 (ongitudinal analysis) 也 叫 纵向 研究 、 纵 贯 面 研究 ， 是 在 比较 长 的 时 间 内 ， 对 一 
个 个 体 或 一 些 个 体 的 心理 发 展 进行 有 系统 的 定期 的 研究 ， 在 心理 学 、 教 育 学 、 社 会 学 、 医 学 
等 学 科 中 被 广泛 使 用 。 随 着 当代 统计 技术 的 不 断 发 展 ， 特 别 是 社会 科学 统计 方法 的 发 展 ， 
用 结构 方程 模型 、 多 层 线 性 模型 等 统计 方法 对 追踪 数据 进行 研究 ， 能 关注 更 多 的 信息 ， 比 如 
增长 趋势 的 个 体 差 异 ， 以 及 影响 个 体 差异 的 各 层级 水 平 的 因素 等 。 由 此 可 见 ， 追 踪 研 究 越 来 
越 细 化 ， 所 关注 的 问题 也 越 来 越 深入 。 

结合 以 往 的 研究 (Usami et al., 2019; Х1% 25, 17%, 2003; SE, FEF, 2019), AY PLS, 
结 出 妃 踪 研究 主要 关注 三 个 方面 的 问题 : (1) 均 值 差 异 的 比较 , (2) 对 多 变量 相互 影响 进行 解释 
以 及 (3) 描 述 总 体 增长 趋势 。 对 于 均值 差异 的 比较 ， 主 要 考虑 到 消除 系统 偏差 、 减 小 误差 方差 
等 因素 ， 可 以 用 重复 测量 的 方差 分 析 或 多 元 方差 分 析 等 方法 。 第 二 类 研究 旨 在 探讨 多 变量 之 
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间 的 往复 式 影 响 (reciprocal relations， 或 相互 影响 )， 与 更 严格 的 因果 影响 (causal effect) 进 行 区 
AY Gi AER, 2017)。 可 以 通过 定义 交叉 清 后 模型 (Cross-Lagged Model, CLM), 探讨 多 变量 变化 
的 相互 影响 。 最 后 一 类 研究 关注 的 是 更 为 宏观 的 增长 趋势 或 增长 趋势 , 引进 潜 增 长 模型 Latent 


Growth Model, LGM). 

在 长 期 的 应 用 当中 ,对 上 述 几 个 问题 的 研究 往往 都 是 孤立 的 。 随 着 近期 一 些 模 型 的 产 出 ， 
研究 者 可 以 按照 研究 问题 ， 同 时 或 选择 性 地 考察 往复 式 影 响 和 增长 趋势 。 从 应 用 者 的 角度 ， 
整合 此 类 的 追踪 研究 将 要 解决 的 问题 ， 给 研究 者 提出 新 的 研究 框架 和 研究 建议 。 


1 模型 定义 


本 节 首 先 介 绍 往复 式 影响 与 增长 趋势 的 两 个 “基础 模型 ”, 随后 引入 变异 分 解 的 具体 议题 ， 


最 后 介绍 结合 了 两 个 基础 模型 的 整合 模型 。 


1.1 交叉 滞后 模型 : 多 变量 往复 式 影响 


探讨 变量 之 间 的 往复 式 影响 常用 的 模型 是 交叉 消 后 模型 (Cross-Lagged Model, CLM), X 


称 为 因果 模型 、 交 又 滞后 面板 模型 (Cross-Lagged Panel Model). 、 自 回归 交叉 滞后 模型 


(Autoregressive Cross-Lagged Model) (Кеппу & Harackiewicz, 1979); 后 来 随 着 结构 方程 模型 
(Structural Equation Modeling, SEM) 技 术 的 发 展 , 能 够 很 容易 地 被 研究 者 所 应 用 (例如 Jöreskog, 


1970)。 交 叉 滞后 模型 以 自 回 归 为 基础 (公式 1), 分 析 所 测 变 量 的 跨 时 稳定 性 (emporal stability, 


参见 Hamaker et al., 2015). 


Хи = Hg + Xa р 


Др, x = В +dy, 。 自 回归 模型 中 ， 观 测 指标 a BRPO i 在 第 1 次 测量 上 的 


fh, АЖ, 为 自 回归 系数 ， 描 述 了 变量 本 身 随 时 间 变 化 的 效应 大 小 ， 即 跨 时 稳定 性 ， 数 值 越 
大 则 稳定 性 越 高 ;dst 为 残 差 ，px 为 方程 的 截 距 。 容 易 拓展 到 多 变量 情形 ， 例 如 含有 两 个 变 
量 的 情况 (公式 2， 图 1)。 


Xi = My + Xir 


Ун Hy, Yir (2) 


ЖН, х, = а-а У Ук 十 dui yi = „улу + У.Х) + а 。 除 了 自 回归 系 


Bl Ba FU В, Д, Ya 和 уу ABC RRL KREE x( 或 变量 y) 在 第 上 次 的 测量 受到 第 1-1 
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次 变量 y( 或 变量 的 影响 ; du 和 а, 表示 残 差 , 服从 联合 正 态 分 布 。 有 研究 者 (Zyphur Allison, 


et al., 2020) 认 为 交叉 滞后 影响 不 应 仅 存 在 于 变量 本 身 ， 


环境 (变量 之 外 的 偶然 因素 ) 也 可 能 对 个 


体 的 后 续 表现 产生 影响 。 故 自 回归 参数 与 交叉 滞后 参数 拓展 到 第 1-1 次 残 差 对 第 1 次 测量 的 影 


啊 ， 被 称 之 为 广义 交叉 滞后 模型 (General Cross-Lagged 


Model; ЖІ 等 , 2021)。 区 别 于 变量 的 


跨 时 稳定 性 ， 残 差 对 下 一 次 测量 的 影响 包含 了 “偶然 因素 ”的 贡献 ， 使 得 模型 能 同时 处 理 变量 
本 身 的 前 后 变化 和 外 在 偶然 因素 对 其 变化 的 影响 。 对 于 传统 CLM 模型 一 般 满足 2 次 的 测量 


就 能 达到 识别 要 求 (Usami et al,，2019);， 而 广义 交叉 滞后 模型 则 需要 至 少 4 次 测量 才能 识别 


(Zyphur, Allison, et al., 2020). 


d, 


研究 者 对 СІМ 模型 中 最 感 兴趣 的 模型 参数 为 自 回 归 系 数 8 RUE SC Js E ІН ЖА y OC 


等 ，201$) 。 此 类 模型 揭示 的 变量 之 间 的 关系 通常 被 研究 者 定义 为 “往复 式 影响 (reciprocal 
effects)”， 即 随 着 时 间 的 推移 ， 每 个 变量 的 变化 都 会 影响 其 他 变量 的 变化 (Usami et al., 2019). 
但 值得 注意 的 是 ， 往 复式 影响 和 因果 影响 仍然 不 能 划 等 号 ， 后 者 需要 满足 更 多 的 前 提 假 设 。 


一 般 地 ,在 严格 的 实验 设计 下 ， 如 果 已 经 控制 了 无 关 3 


Sit, CLM 可 以 用 来 揭示 变量 之 间 的 


TE AG ER, 2017). 


CLM 模型 定义 相对 单 , 对 测量 次 数 要 求 也 不 多 , 也 能 在 一 定 程度 上 揭示 不 同 变量 随时 间 


变化 的 相互 影响 关系 , 应 用 领域 涉及 到 个 体 特质 教育 成 果 与 精神 疾病 特征 等 领域 的 长 期 性 、 


滞后 性 影响 的 追踪 研究 (Cacioppo et al., 2017; Chiu & Du, 2019; Liu & Hou, 2018; Selzler et al., 


2019)。 国 内 这 方面 研究 主要 集中 在 青少年 问题 行为 、 


AE 等 ,2019; “т 等 , 2018)。 


抑郁 、 能 力 发 展 的 影响 因素 等 方面 ( 陈 
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1.2 潜 增 长 模型 : 个 体 增长 趋势 


另 一 类 问题 是 描述 个 体 增 长 趋势 , 包含 总 体 增长 形态 的 特点 及 其 是 否 存在 个 体 差异 ( 即 发 
展 群 体 的 同 质 性 )。 潜 增长 模型 (Latent Growth Model, LGM, Bollen & Curran, 2006) 是 建立 在 
SEM 视角 下 ,通过 定义 带 有 均值 结构 的 验证 性 因素 分 析 模 型 ， 固定 特殊 载荷 大 小 来 实现 增长 
趋势 的 描述 。 以 线性 LGM 为 例 (公式 3- 公 式 4): 


测量 部 分 六 =m + @ —Drp; +, (3) 


UO 

其 中 ，yi 表 示 被 试 i TESS 上 次 的 观测 指标 ， 其 测量 误差 分 别 对 应 sa ЖУА TERES 
ЖӘН qu 和 线性 斜率 mi 潜 变 量 因子 的 均值 分 别 对 应 кү Жік» 方差 对 应 于 va) var), Gi 
和 ;服从 联合 正 态 分 布 。 其 中 , 潜 变 量 m RU qo 的 系数 矩阵 A 可 以 根据 时 间 和 形态 灵活 定义 
或 增 减 斜 率 因 子 ， 诸 如 二 次 增长 模型 、 不 定义 增长 曲线 模型 和 单 因 素 增 长 曲线 模型 等 ( 唐 文清 

等 , 2014; 王 备 成 , 毕 向 阳 , 2018)。 除 了 连续 的 增长 之 外 ， 还 可 以 定义 间断 增长 的 趋势 ， 如 多 
阶段 增长 模型 (piecewise growth model, Liu et al., 2018; Lock et al, 2018)。 一 般 对 于 两 个 因子 的 
线性 增长 模型 ，3 次 测量 是 保证 模型 识别 的 最 低 要 求 ， 而 如 果 需 要 增加 增长 因子 (如 二 次 、 多 
阶段 因子 )， 需 要 增加 测量 次 数 。 


图 2 潜 增长 模型 示意 图 
描述 个 体 增长 趋势 有 助 于 研究 者 把 握 个 体 发 展 的 形态 ， 从 而 根据 不 同 阶段 的 发 展 特 征 制 


定 不 同 的 教育 政策 或 干预 手段 。 如 果 需 要 从 政策 或 根源 上 来 干预 增长 ， 模 型 中 一 般 会 考虑 非 
时 变 (time-invarianb 协 变量 ， 他 们 不 会 随 着 时 间 的 变化 而 变化 ， 其 影响 也 会 较为 持久 、 沸 后 。 
如 学 生 的 性 别 、 社 会 经 济 地 位 、 历 史 事 件 等 相对 稳定 的 因素 (Diallo et al., 2017; Liu et al., 2020). 


4 


是 及 时 的 、 变 化 的 。 此 时 ， 多 水 平 模型 (Multilevel Model，MILM) 中 的 时 变 协 变量 定义 方式 也 


可 以 纳入 LGM 的 研究 当中 (Liu etal., 2016; МсСоасһ & Kaniskan, 2010). 

LGM 模型 定义 也 相对 简单 ， 研 究 者 如 果 有 3 次 以 上 的 追踪 则 可 使 用 该 模型 建 模 。LGM 
将 变化 分 解 为 初始 状态 和 增长 趋势 , 可 以 宏观 地 从 两 个 维度 来 考察 总 体 增长 趋势 , 在 SEM HE 
架 下 定义 也 相当 灵活 。 国 内 的 应 用 研究 近年 逐日 增多 , 集中 在 儿童 青少年 社会 能 力 、 孤独 感 、 
自我 概念 、 积 极 适 应 的 发 展 和 老年 人 的 认 知 功能 发 展 等 方面 ( 侯 桂 云 等 , 2018; 李 彩 娜 等 ， 
2017; 李 振 兴 等 , 2017; 刘 俊 升 等 , 2013; 张 晓 , 2011)。 其 中 ， 也 包含 时 间 跨 度 较 短 的 状态 类 
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变量 ， 如 人 际 冒 犯 、 抑 郁 等 ( 程 刚 等 , 2016; 杨 满 云 55,2017). 


1.3 测量 误差 与 随机 截 距 


以 上 回顾 了 两 个 “基础 模型 "可 以 在 他 们 基础 之 上 添加 一 些 特殊 议题 ， 将 变异 进一步 分 
解 。 以 交叉 滞后 模型 为 载体 ， 借 用 SEM 的 思路 ， 考 察 每 个 指标 的 测量 误差 ， 可 以 建立 “因子 


交叉 滞后 模型 (Factor Cross-Lagged Model, FCLM)"(Jóreskog & Sórbom, 1974; McArdle, 2009; 


Usami et al., 2015). 


it xit xit 
Уи = Fs 十 € ji; (5) 
c 结构 部 分 : 


Fu = Ha + P ed) + Vaf jian ға, 
К = Hy + „ЕЁ л) ЕЗ yy ааа) + dy, (6) 


Ж, Fa 和 Ру 表示 潜 变 量 ， 是 分 别 排除 了 测量 误差 ew 和 si 之 后 的 真 分 数 。 在 ЕСІМ 


模型 中 ， 交 叉 滞 后 系数 都 建立 在 潜 变 量 上 ， 当 测量 误差 为 0 时，FCLM 可 以 简化 成 标准 CLM 


模型 。 与 此 同时 ，FCLM 当中 每 个 测量 指标 到 因子 之 间 的 路 径 系 数 限定 为 1， 测 量 误差 服从 


ІШ 
ІШ 


联合 正 态 分 布 ， 需 要 至 少 3 次 测量 能 满足 模型 识别 的 要 求 。 如 果 每 次 的 测量 模型 为 多 指标 ， 


可 以 按照 SEM 模型 识别 的 要 求 建立 测量 模型 。 
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图 3 因子 交叉 滞后 模型 示意 图 


另 一 方面 ， 自 回归 系数 虽然 描述 了 跨 时 稳定 性 ， 但 是 其 即使 在 很 高 的 情况 下 ， 经 过 足够 
长 的 时 间 ， 原 始 的 秩序 也 会 丢失 (如 图 1 中 加 > 总 的 自 回 归 间 接 效 应 为 乘积 A opa na, ЖИН 
小 于 1)。 随 机 截 距 交 又 灌 后 模型 (Random Intercept Cross Lagged Model，RI-CLM) 的 建 模 方式 
是 将 交叉 滞后 模型 增加 一 个 随机 截 距 项 ， 用 来 描述 真正 的 非 时 变 的 特质 (tme-invariant 


trait-like stability， 参 见 Hamaker etal., 2015; Mulder & Hamaker, 2020)。 具 体 建 模 方式 如 下 


ж 
Хи = H, Ты T Xi 


Yir — Hy ty Xie 


* * * * _ ж ж уде ZX RES Ж) 
其 中 x, = бл, ЖУУ) * Ga Ya = Bu Vien уа) + d, o 在 这 个 模型 中 ， 


т 和 myi; 表 示 非 时 变 稳定 因子 ， 用 以 呈现 跟 特 质 相 关 的 因素 ,固定 其 对 于 所 有 测量 指标 的 影 
响 都 为 1。 这 里 ， 区 别 于 LGM BUS, gu 和 iy; 均 值 为 0， 其 方差 协 方差 矩阵 Ф 和 截 距 项 yw 
All uy, 自由 估计 。 一 般 要 求 至 少 3 次 测量 能 满足 模型 识别 的 要 求 。 

由 于 引入 了 特质 变量 ， 可 以 将 个 体 发 展 的 变异 分 解 为 个 体 间 变 异 (between-person 
difference) 和 个 体内 变异 (within-person difference)， 这 是 此 类 模型 进行 拓展 的 重要 出 发 点 
(Bainter & Howard, 2016; Curran & Bauer, 2011; Hamaker et al., 2015; Voelkle et al, 2014; Ж 
等 , 202D)。 其 中 ， 个 体 间 变异 来 源 为 提取 所 有 个 体 在 数 次 测量 之 间 的 公共 因素 ， 表 示 个 体 跨 
时 间 的 稳定 的 特征 ;个 体内 变异 是 每 个 个 体 在 数 次 测量 之 间 的 滞留 影响 ， 表 示 个 体 本 次 状态 
受到 先前 状态 的 影响 (Hamaker et al; 2015). ТЕКІ-СІМ 中 ， 加 入 随机 截 距 的 目的 在 于 将 稳定 
的 、 不 随时 间 变 化 的 变异 真正 分 解 为 非 时 变 的 特质 (个 体 间 变异 ) 和 跨 时 稳定 性 (个 体内 变异 ) 
因素 。 此 时 , 自 回归 系数 ps 和 By, 被 称 为 个 体内 滞留 参数 (within-person carry-over, Hamaker et 
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al., 2015), 亦 被 称 为 个 体 惰性 (inertia, Mulder & Hamaker, 2020)。 如 果 将 RLCLM 的 gi ЖІ my; 
限定 成 0 时 ， 就 成 为 一 个 标准 СІМ 模型 ， 表 示 没 有 跨 时 间 的 稳定 特质 ， 即 不 存在 个 体 间 差 


图 5. 随机 截 距 交叉 滞后 模型 示意 图 


在 SEM 框架 下 ， 很 容易 将 测量 误差 与 随机 截 距 模 型 结合 ， 得 到 “特质 -状态 -误差 ”模型 
(Trait-State-Error, TSE, Kenny & Zautra, 1995, 200D)。 并 且 ， 这 个 考察 单个 变量 变化 的 模型 
很 快 拓展 到 多 元 情形 (Luhmann et al., 2011; Zautra et al, 1995)。 首 先 ,TSE 模型 考虑 测量 误差 ， 
即 每 个 题目 都 有 各 自 的 测量 信 度 (公式 5)。 分 解 出 测量 误差 sa 和 syi 后 ， 真 分 数 ҒаЖ Ру BJ 
变异 继续 被 分 解 成 两 个 部 分 :一 部 分 是 稳定 特质 (yi ЖІ ту) ЕРЕ, 一 部 分 是 状态 (Fxa FU F) 
之 间 的 关系 (包括 自 回 归 系 数 pB M В. ELME RÉ ya 与 yy 两 部 分 ): 

Fa = Hha +. + BaPa- + Variety + dus 


Fi = Hy thy T yt Fee 1) + yy Fac) +d 


(8) 
一 般 地 ， 特 质 与 状态 的 影响 都 固定 为 1。 如 果 假定 时 变 的 参数 ， 比 如 自 回归 参数 p、 交 叉 
滞后 参数 》 残 差 和 协 方差 矩阵 都 随 着 c 的 不 同 而 自由 估计 ， 要 求 至 少 4 次 测量 能 满足 模型 | 
别 的 要 求 。 如 果 假 定 非 时 变 的 参数 ， 即 上 述 参数 都 限定 为 相等 ，3 次 测量 即 可 满足 模型 识别 
要 求 。 


LI 


“ 


` 
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图 6 特质 -状态 -误差 模型 示意 图 


TSE 中 ,随机 截 距 的 ты ЯП 7yi 代 表 非 时 变 因子 ， 即 个 体 的 特质 (traib， 是 “个 体 间 差异 ” 


自 回 归 和 交叉 滞后 部 分 表示 了 当前 的 测量 受到 前 一 次 测量 的 影响 ， 是 “个 体内 差异 ”， 此 二 者 


5j КІСІМ 相同 。 与 此 同时 ， 考 察 每 个 测量 指标 的 真 分 数 ， 
被 解释 的 部 分 是 误差 (error)， 此 与 FCLM 相同 。 不 难看 出 ， 


被 称 之 为 个 体 的 状态 (state)， 不 能 
TSE 同时 包含 了 随机 截 距 和 测量 


误差 ,是 将 RECLM 和 FCLM 进行 了 整合 。 也 有 研究 者 将 误差 固定 为 0 后 将 其 看 成 是 考虑 了 


测量 信和 度 的 RI-CLM 模型 (Mulder & Hamaker, 2020). 
该 模型 在 后 续 发 展 中 改名 为 “稳定 特质 自 回归 特质 与 状 


Trait and State Model, STARTS)” (Kenny & Zautra, 2001; Schmitt & Steyer, 1993)。STARTS 在 


态 模型 (Stable Trait Autoregressive 


命名 上 区 分 了 两 种 “特质 ”: 稳定 特质 指 来 自 于 随机 截 距 的 部 分 Ca 和 туі)» МАНІ ЫР; 
自 回归 特质 (ps 和 0 即 “ 个 体内 滞留 影响 ”。 状 态 指 每 次 测量 的 真 分 数 (Esa 和 Ға), BAER M 
量 误差 >。 所 以 在 命名 上 ，STARTS 更 为 精确 。 由 于 STARTS 和 TSE 只 是 命名 上 的 区 分 ， 为 


了 简化 ， 仍 然 将 此 模型 称 之 为 特质 -状态 -误差 ”模型 。 


TSE 模型 优点 在 于 可 以 分 离 出 状态 和 特质 两 个 因子 的 心理 学 概念 ， 同 时 ， 它 将 往复 式 影 
响 精 确定 位 到 个 体内 层面 ， 被 广泛 地 应 用 在 发 展 心理 学 和 临床 领域 中 。 如 抑郁 ， 可 分 为 特质 
与 状态 两 个 层面， 前 者 是 个 体 稳定 的 特质 ， 随 时 间 推移 保持 稳定 ， 而 后 者 是 情境 变量 ， 随 时 
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间 推 移 而 变化 (Hazel & Hankin, 2014; Kiken et al., 2015; Masselink et al., 2018)。 分 离 出 稳定 的 
特质 和 变化 的 状态 之 后 ， 研 究 者 可 以 根据 其 占 比 大 小 ， 进 行 适当 的 干预 ， 特 别 是 从 特质 因素 
的 根源 出 发 ， 寻 找 更 有 效 的 诊断 方法 。 例 如 IMasselink 等 人 (2018) 通 过 三 个 不 同样 本 、 不 同时 
间 跨 度 的 研究 发 现 ，RLECLM 建 模 指向 “个 体内 仅 有 自尊 负 向 预测 抑郁 ”的 结论 ， 自 尊 和 抑郁 
的 负 相 关 是 源 于 特质 层面 ， 而 CLM 建 模 则 指向 “个 体内 自尊 和 抑郁 之 间 存 在 负 向 的 往复 式 影 
响 ” 的 结论 。 该 研究 的 差异 化 结果 恰好 证 实 了 贝克 认 知 理论 (Beck’s cognitive theory) 当 中 的 易 
感 模 型 (vulnerability model): 低 自 苯 的 个 体 容易 抑郁 ， 而 本 来 就 难以 检测 的 创伤 模型 (scar 
modeD) 在 RI-CLM 模型 中 被 证 伪 。 说明 使 用 CLM 检测 的 灵敏 度 不 如 RI-CLM， 后 者 能 更 精准 
地 定位 出 差异 源 于 何 处 。 


1.4 自 回 归 模 型 与 潜 增 长 模型 结合 


将 交叉 滞后 模型 与 潜 增长 模型 结合 ， 可 以 同时 考察 往复 式 影响 与 增长 趋势 ， 其 建 模 方式 
主要 有 两 种 思路 。 第 一 种 建 模 方法 是 以 交叉 滞后 模型 为 基础 模型 ， 在 其 基础 之 上 增加 潜 增 长 


IF. 其 中 较 早 提出 来 的 是 Curran 和 Bollen (2001) 的 自 回 归 潜 增长 模型 (Autoregressive Latent 


Trajectory, ALT, Curran & Bollen, 2001)。 其 建 模 方式 如 下 。 

X, =Ma +0 0р Bax + Y yi + dat D 

y — hy t -Dhat E қат ауаға, 20 (9) 

在 结合 的 时 候 , 注意 到 x ЛП у, 定义 的 区 别 。 AVES i RS, хі 和 yi 是 不 受 模 
型 内 其 他 变量 影响 的 ， 属 于 “外 源 变量 *"， 而 在 定义 增长 因子 的 时 候 ，xi 和 六 被 当成 “指标 ”。 
因为 ALT 的 载体 为 交叉 滞后 模型 ， 所 以 需要 将 立 和 考虑 成 外 源 变量 ， 即 设 定 立 和 yi 与 增 


长 因子 之 间 为 相关 并 自由 估计 (Curran & Bollen, 2001; Jongerling & Hamaker, 2011; Zyphur, 


Voelkle et al., 2020)。 此 外 ,如 果 将 斜率 项 ур: 和 myi Е 0, 则 得 到 RI-CLM 模型 (公式 7), 
故 RI-CLM Ж ALT 的 一 个 特例 。 该 模型 需要 有 5 次 测量 才能 满足 时 变 参数 模型 的 识别 要 求 ， 
而 如 果 限 定 为 非 时 变 参 数 ，4 次 测量 即 可 。Usami 等 人 (2019) 也 把 ALT 模型 称 之 为 “累积 增长 
模型 ”， 增 长 因子 被 看 成 是 两 个 “累积 因子 (accumulatingfactors)”。 
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图 7 自 回 归 潜 增长 模型 示意 图 


第 二 种 建 模 方式 以 潜 增 长 模型 为 基础 模型 ， 在 其 基础 之 上 增添 交叉 滞后 的 影响 系数 。 其 


代表 为 结构 化 残 差 潜 增长 模型 (Latent Curve Model with Structured Residuals, LCM-SR, Bainter 


& Howard, 2016; Curran et al., 2014). E49 ALT 模型 的 出 发 点 不 同 : ALT 出 发 点 是 自 回 归 ， 


以 此 添加 个 体 间 变 异 ( 从 个 体内 影响 中 分 离 出 个 体 间 变 异 ); LCM-SR 出 发 点 是 潜 
添加 交叉 滞后 影响 (从 个 体 间 变 异种 分 离 出 个 体内 影响 )。 所 以 ，LCM-SR 不 纠结 


首长 ， 以 此 


第 一 次 测 


量 是 否 为 外 源 变 量 ， 直 接 将 因子 潜 变 量 作用 于 所 有 测量 指标 ， 并 且 在 测量 部 分 添加 交叉 滞后 


效应 (图 8)。 建 模 可 以 表示 为 
X, = Tha + E-D + Жағы 十 ПЕТ td, 


Vie = Ny + (t — Ty, + муа) T Уу шй р dy, (10) 


RP, iel sicD 是 不 能 被 增长 因子 所 解释 的 误差 项 , 在 他 们 之 间 建 立 自 回归 和 交叉 
滞后 参数 。 在 LCM-SR 中 ,更 清晰 地 看 到 观测 变量 的 变异 被 分 解 到 了 被 试 间 和 被 试 内 两 部 分 。 


被 试 间 变异 表示 为 潜 增 长 模型 部 分 ， 被 截 距 项 (74s、wyi) 和 斜率 项 (1pwx、112yi) 解 释 ; 


被 试 内 变 


异 是 误差 项 grow 和 aie， 将 自 回归 和 交叉 灌 后 的 影响 直接 定义 在 测量 误差 之 间 (Bainter & 
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Howard, 2016; Mulder & Hamaker, 2020) 。 由 于 误差 项 本 身 也 是 潜 变 量 ， 故 LCMLSR 模型 的 组 


内 差异 建立 在 潜 变 量 上 。 
实际 中 如 何 定义 误差 之 间 的 关系 ? 如 果 将 e* 看 成 是 FY 的 一 个 特例 , 那 e* 即 潜 变 量 因子 ， 


定义 其 观测 指标 x 和 y 的 “测量 误差 ?为 0( 观 测 指 标 不 再 含有 其 他 的 误差 因子 )， 就 巧妙 地 将 理 


论 定 义 操作 化 了 (参见 Mulder & Hamaker 2020; Zyphur Allison, et al., 2020)。LCM-SR 要 求 4 


次 测量 满足 时 变 参 数 的 识别 要 求 ， 如 果 限 定 为 非 时 变 参数 则 需要 З 次 测量 


对 比 两 种 建 模 思 路 可 知 ,， 除了 是 否 假设 第 一 次 测量 为 外 源 变量 之 外 ，ALT 和 LCM-SR 的 
建 模 差 异 集中 在 交叉 滞后 的 影响 是 哪 一 类 变量 所 指向 的 。ALT 当中 的 交叉 滞后 影响 直接 建立 


在 观测 变量 上 ， 即 不 考虑 测量 误差 ;而 LCM-SR 当中 的 交叉 滞后 影响 建立 在 测量 误差 上 ， 


潜 变 量 上 。 如 果 将 ALT 模型 进行 拓展 ， 考 虑 其 每 个 观测 变量 的 测量 模型 ， 但 是 限定 这 些 测 量 


р 


E. 


模型 的 “测量 误差 ?为 0, 即 可 得 到 LCM-SR 模型 ( 除 第 一 次 测量 为 外 源 变量 )。 换 言 之 , LCM-SR 


可 以 粗略 看 成 是 一 个 带 有 因子 的 ALT 模型 在 经 过 限定 其 测量 误差 为 0 的 基础 上 得 到 的 。 


8 结构 化 残 差 潜 增 长 模型 示意 图 
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除 上 述 两 类 建 模 思想 以 外 ， 潜 变化 分 数 模 型 (Latent Change Score Model, LCS), WIRA 


异 分 数 模型 Latent Difference Score Model，LDS)， 也 是 同一 时 期 衍生 出 来 的 同时 关注 交叉 漳 


后 与 潜 增 长 的 模型 (McArdle & Hamagami, 2001). LCS 在 模型 定义 上 , 首先 纳入 了 测量 误差 ( 公 
式 5)， 然 后 构造 “因素 差异 "(公式 11), 并 且 用 差异 分 数 作为 指标 构造 增长 特质 (ps 和 ту» A 
式 12)， 即 

Fa = Fay + АЁ 


xit 
F Ёш) +АЕ (11) 


yit 一 
AF. = AN xi TIE d 十 KP nan 本 d ut e D 


AF i = ay yi FRE efi Ед Vo age ға, > 1) 


(12) 

模型 中 的 增长 特质 my 和 yy; 被 称 为 “斜率 "， 但 容易 将 上 述 公 式 进行 转换 ， 将 公式 12 带 
入 公式 11， 得 到 

Fa = 0. ++ Ba) + Уы Py +d,,(t >D 


Fi, = 0; 十 (1+ Pb) Erian 23 P Eun 十 ам >) 


(13) 
可 知 , LESCAR 13) 中 的 “斜率 因子 一 一 Nw 和 jy 一 一 其 实 本 质 上 仍然 是 ALT 当中 的 一 


Y 个 “ 截 距 因子 ”(Usami et al., 2015; Usami et al., 2016)。 因 子 载荷 wk 和 ow 代表 了 每 次 测量 的 权 
© 重 。 在 传统 LCS ANP, о, 和 оу 是 设 定 为 非 时 变 的 ， 即 wu=oy=1; 若 将 此 限定 放宽 ， 即 得 
到 时 变 的 截 距 项 ошл ЖІ oy iyi 也 被 称 为 三 重 变化 分 数 模型 (triple change score, TCS, McArdle 
& Nesselroade，2014)。 这 样 就 可 以 在 LCS 基础 上 对 协 方 差 的 进行 限定 ， 得 到 ЕСІМ 模型 ( 公 
式 6)。 所 以 ,二 者 具备 从 套 关系 ,可 将 FCLM 看 成 是 LCS 的 一 个 特例 (Usami et al., 2015, 2016). 
除 此 之 外 ， 二 者 的 自 回 归 系 数 的 定义 也 略 有 差异 。 

同时 关注 潜 增 长 和 交叉 灌 后 的 模型 在 模型 的 解释 与 应 用 上 相对 复杂 ， 常 被 应 用 在 发 展 与 
临床 研究 当中 (Berry & Willoughby, 2017; Cole et al., 2005; Curran et al., 2014; Ding et al., 2020; 
Malone et al., 2004)。 这 类 模型 解决 的 核心 问题 是 ， 控制 了 稳定 的 个 体 间 变 异 以 后 ， 个 体内 是 
否 仍然 存在 发 展 惰性 或 变量 之 间 的 往复 式 影响 。 由 于 在 СІМ 中 往复 式 影响 杂 焰 了 稳定 的 成 
分 ， 所 以 传统 模型 的 参数 估计 与 整合 后 的 模型 结果 或 大 相 径 庭 ， 而 如 果 考 察 的 目标 变量 确实 
存在 这 种 “特质 (研究 中 常 提 到 的 抑郁 、 焦 虑 、 自 尊 等 心理 学 成 分 )， 不 经 控制 的 模型 会 高 估 
这 种 影响 。 例 如 , Berry 和 Willoughby(2017) 的 研究 中 通过 两 组 真实 数据 对 这 类 数据 进行 建 模 ， 
探索 儿童 的 攻击 行为 与 父母 的 体罚 ( 打 屁股 ，spanking) 之 间 的 关系 。 发 现在 纳入 了 增长 模型 之 
后 ， 体 罚 行 为 随 着 时 间 推 移 ， 频 次 降低 (斜率 为 负 )， 但 是 攻击 行为 有 所 增加 ;两 种 行为 的 交 


又 滞后 影响 不 显著 。 但 是 对 比 传统 CLM,， 两 种 行为 之 间 存 在 交叉 滞后 影响 。 研 究 者 指出 ， 
12 


B 


复式 影响 或 不 应 发 生 在 每 个 层面 。 因 为 个 体 间 的 变异 应 是 稳定 的 非 时 变 效应 ， 而 个 体内 变异 


应 是 变化 的 、 波 动 的 ， 所 以 剔除 个 体 间 变异 对 于 模型 的 可 靠 性 更 是 必要 的 前 提 (Berry & 


Willoughby, 2017)。 所 以 ， 在 个 体 间 层 面 ， 儿 童 的 攻击 行为 受到 父母 体罚 行为 的 正 向 影响 ， 然 


而 这 并 不 会 发 生 在 个 体内 。 


2 模型 整合 与 拓展 


这 类 研究 的 核心 问题 在 于 如 何 分 解 变异 。 潜 增长 模型 所 提取 出 来 的 总 体 增长 趋势 被 看 成 


2016; Hamaker et al., 2015; Murayama et al., 2017)。 在 提取 公共 因子 之 前 ， 往 复式 影响 发 生 在 


是 个 体 间 变异 ;不 能 被 它 解释 的 ， 即 交叉 滞后 模型 探索 的 ， 是 个 体内 变异 (Bainter & Howard, 


个 体 的 所 有 层面 ， 但 是 当 分 离 出 特质 变量 之 后 ， 往 复式 影响 被 定位 到 个 体内 层面 。 本 节 的 模 


型 整合 与 拓展 思路 将 按照 追踪 研究 主要 解决 的 问题 来 展开 ， 
响 和 增长 趋势 的 问题 。 


则 可 建立 潜 增 长 模型 (LGM)。 然 后 ， 分 别 以 CLM 3E LGM 为 基础 模型 ， 添 加 相应 的 变异 分 解 


旨 在 解决 如 何 精确 定位 往复 式 影 


首先 , 研究 者 若 关注 往复 式 影响 , 可 以 建立 交叉 滞后 模型 (CLMD); 若 关 注 个 体 增长 趋势， 


成 分 。 诸 如 测量 误差 ， 以 潜 变 量 指标 代替 观测 变量 进行 建 模 ， 或 者 根据 研究 需要 只 关注 截 距 


因子 或 斜率 因子 。 若 考察 所 有 因素 ， 可 以 将 上 述 模型 进行 拓展 ， 得 到 整合 的 模型 一 一 以 CLM 


为 基础 模型 的 潜 变 量 自 回 归 洪 增长 模型 Latent Variable Autoregressive Latent Trajectory, 


LV-ALT, Bianconcini & Bolen，2018)， 以 及 以 LGM 为 基础 模型 的 因子 结构 化 潜 增 长 模型 


(Factor Latent Curve Model with Structured Reciprocals, FLCM-SR). 
以 LVALT СІМ 为 基础 模型 )， 其 模型 建构 表示 如 公式 14 所 示 ( 图 9). 


Fi =7 tl -Dn +f, F pm у ub. ей TUI 


F 


yit 
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9 潜 变 量 自 回归 潜 增 长 模型 示意 图 


建 模 当 中 的 各 个 参数 同 前 文 。 首先 考虑 测量 误差 (公式 5)， 自 回归 和 增长 曲线 均 在 真 分 数 
Fa 和 Fa 上 建 模 。 其 次 ， 考 虑 个 体 的 增长 形态 ， 即 截 距 Gps 和 туу ЯН (ра 和 қу)» HA 
数 固定 为 特定 的 增长 形态 (其 与 第 一 次 测量 真 分 数 Fa 和 Fy 1 设 定 为 相关 而 非 增长 形态 )。 在 此 
基础 上 ， 建 立 自 回 归 影 响 (B64，py) 和 往复 式 影 响 (y4，yyt)。LV-ALT 模型 识别 要 求 可 分 两 阶段 
BET: 结构 部 分 的 识别 同 ALT， 时 变 效应 需要 5 次 测量 ，4 次 测量 要 求 限 定 线 性 增长 ，3 次 
测量 要 求 进一步 限定 自 回归 系数 相等 ; 测量 部 分 的 识别 与 一 般 验 证 性 因素 分 析 模 型 相同 
(Bianconcini & Bollen, 2018). 

另 一 种 建 模 思路 ， 以 LGM 模型 为 基础 模型 ， 包 含 了 所 有 分 解 变异 成 分 的 模型 即 为 因子 
结构 化 潜 增 长 模型 (FLCM-SR)。 其 构建 方法 如 公式 15. 

Faha t C-ie “Б.Ғаз%7 +, 
F a — yt T Dist „Елу + Veh ign + d 


yit (15) 
这 里 ， 观 测 指 标的 变异 首先 被 看 成 是 包含 真 分 数 和 测量 误差 (公式 5)， 然 后 将 真 分 数 Fyi 
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RIF, UAE SEA MEC AU СӨЗІН RI m i по ЖП oy) MENERA Py 
和 _Fsiti) 的 自 回归 影响 和 往复 式 影响 ) 和 残 差 ls 和 di)。 注 意 到 公式 15 中 的 PATAR 
10 中 的 s*， 这 意味 着 FCLM-SR 模型 中 将 结构 化 残 差 因子 看 成 是 真正 的 潜 变量 ， 追 加 每 个 观 
测 指标 真正 的 测量 误差 (公式 5)。 具 体操 作 层面 ， 在 LCM-SR 基础 之 上 放宽 观测 指标 测量 误 
差 为 0 的 假设 即 可 (附录 3)。 


10 因子 结构 化 潜 增 长 模型 示意 图 
对 比 LV-ALT， 可 知 FCLM-SR 将 公式 14 中 1 的 限定 放宽 ， 忽 略 第 一 次 测量 的 外 源 变 
量 性 质 。 从 图 10 来 看 ， 截 距 mi ту 直接 作用 于 xi 和 y 并 固定 载荷 为 1， 而 不 再 添加 因子 
之 间 的 相关 ; 其 他 的 建 模 与 LV-ALT 相同 。 换 言 之 ， 两 类 模型 在 建 模 过 程 中 的 核心 问题 仍然 
为 “是 否 将 第 一 次 测量 看 成 外 源 变量 *， 其 他 的 参数 设 定 都 只 是 在 整体 框架 当中 的 特殊 案例 。 


这 样 的 假设 在 实证 当中 具有 可 操作 性 ， 可 以 让 第 一 次 观测 指标 从 属于 个 体 间 变异 ， 即 从 研究 
一 开始 就 存在 固定 的 特质 。 因 为 实证 研究 的 首次 测量 点 并 不 一 定 是 特质 发 展 当中 真正 的 “ 原 
点 ”, 或 存在 于 整个 研究 外 部 的 真正 的 外 源 变量 。 故 可 在 实证 过 程 中 直接 用 第 一 次 的 观测 指标 
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建 模 了 ， 由 此 一 来 对 于 结构 模型 识别 的 要 求 也 可 以 减少 一 次 观测 指标 (时 变 模 型 4 次 测量 、 非 
时 变 模 型 3 次 测量 )。 

Usami 等 人 (2019) 在 对 上 述 模 型 进行 整合 的 时 候 ， 将 变异 分 解 为 : (1) 测 量 误差 
(measurement erroD， 诸 如 因子 交叉 滞后 模型 、 特 质 -状态 -误差 模型 中 对 测量 信 度 的 考察 ; Q) 
成 分 方程 (decomposition equation)， 即 潜 增 长 部 分 ， 包括 只 有 随机 截 距 与 同时 包含 截 距 和 斜率 

量 看 成 外 源 变 量 之 后 的 建 模 。 但 是 该 整合 框架 
问题 在 于 ， 成 分 方程 和 动态 方程 并 不 能 够 同时 建 模 一 一 实质 上 ， 动 态 方程 表示 了 以 交叉 滞后 
模型 为 基础 的 建 模 思想 ， 而 成 分 方程 表示 了 以 潜 增 长 模型 的 建 模 思 想 。 故 重新 对 多 变量 追踪 
究 的 系列 模型 所 包含 的 成 分 进行 分 解 ， 分 别 以 交叉 滞后 模型 和 潜 增 长 模型 为 基础 模型 ， 直 
接 对 研究 者 都 比较 熟悉 的 几 个 成 分 进行 增 减 建 模 。 拓 展 框 架 如 表 1 所 示 。 


| 


的 情况 ; (OIA) ЁЕ(Чупапис), 183826 — Uc 


= 


zn 


表 1 交叉 滞后 与 潜 增长 系列 模型 的 成 分 构成 一 览 
ВИН ЖАС REAT WERA MER% 


基础 模型 1: 交叉 滞后 模型 (CLM) ү Е 2 4 : 
因子 交叉 滞后 模型 EC-CLMD) 4 | Xl 4 
ЕЛЕЕ S SC Je ER (RI-CLM) Y a 4 
特质 -状态 -误差 模型 (TSE/STARTS) y RE - N J 
自 回归 潜 增长 模型 (ALT) RE ү» у : 
潜 变 量 自 回 归 潜 增 长 模型 (LV-ALT) ү» ИШ 4 N 


基础 模型 2， 潜 增长 模型 (LGMD = y Р : 
结构 化 残 差 潜 增 长 模型 (LCM-SR) = ү ү ү RE 

因子 结构 化 残 差 潜 增长 模型 (FLCM-SR) RI ү \ ү 

FE 截 距 因子 的 均值 为 0。"” 假 设 第 一 次 观测 变量 为 外 源 变量 (与 增长 因子 为 相关 )。“ 限定 为 0。 


交叉 滞后 模型 当中 ， 不 考虑 测量 误差 ， 也 不 考虑 增长 形态 ， 建 模 最 简单 。 在 此 基础 上 ， 
可 以 分 别 考察 测量 误差 (因子 交叉 滞后 模型 ) 和 随机 截 距 (随机 和 截 距 交 叉 滞 后 模型 ), 或 将 二 者 同 
时 考虑 (特质 -状态 -误差 模型 )， 在 此 基础 上 ， 加 入 增长 模型 ( 自 回归 潜 增 长 模型 ) 及 其 测量 误差 
( 潜 变 量 自 回 归 潜 增长 模型 )。 另 一 个 思路 是 潜 增 长 模型 ， 不 考虑 交叉 滞后 ， 也 不 考虑 测量 误 
差 ， 建 模 也 相对 简单 。 在 此 基础 上 ， 考 虑 交叉 滞后 成 分 (结构 化 潜 增 长 模型 )， 或 进一步 考虑 
测量 误差 (因子 结构 化 残 差 潜 增 长 模型 )。 两 种 思路 最 终 统一 于 LV-ALT/FCLM-SR. WREE 
选用 模型 的 时 候 ， 应 该 以 某 一 类 模型 为 基础 模型 ， 在 此 基础 上 添加 相应 的 变异 分 解 成 分 。 
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3 实证 应 用 : 往复 式 影 响 与 个 体 发 展 的 模型 比较 


采用 美国 国家 教育 统计 中 心 (National Center for Education Statistics，NCES) 提 供 的 早期 儿 


童 的 追踪 研究 -幼儿 园 版 (Early Childhood Longitudinal Study-Kindergarten cohort, ECLS-K, 


https:/nces.ed.gov/ecls/)。 该 追踪 研究 调查 了 1998-1999 年 入 学 的 21409 名 学 生 ， 用 长 达 9 年 
的 追踪 调查 了 从 幼儿 园 到 八 年 级 的 学 生 的 智力 发 展 情况 和 来 自 不 同 层面 的 影响 因素 。 数 据 测 
量 了 7 个 时 间 点 : 幼儿 园 秋 季 学 期 ， 幼 儿 园 春季 学 期 ， 一 年 级 秋季 学 期 ， 一 年 级 春季 学 期 ， 
三 年 级 春季 学 期 ， 五 年 级 春季 学 期 和 八 年 级 春季 学 期 。 测 量 的 指标 有 直接 认 知 能 力 (包括 阅读 
成 绩 、 数 学 成 绩 和 科学 成 绩 )、 间 接 认 知 能 力 (包括 教师 对 学 生 认 知 能 力 的 整体 评价 ， 采 用 行 
量 表 计 分 )、 学 生 的 自我 描述 以 及 学 生 的 自我 概念 和 控制 点 。 由 于 第 三 次 (一 年 级 下 期 ) 只 对 
其 中 30% 的 被 试 进行 了 调查 ， 本 研究 中 没有 采用 这 一 次 的 调查 数据 ， 只 使 用 了 其 余 六 次 测量 
的 数据 ， 有 效 样本 21049 名 。 第 一 次 测试 时 学 生 的 平均 年 龄 68.47 月 ， 标 准 差 4.21 月 ， 最 大 
值 79 个 月 ， 最 小 值 54 个 月 。 测 量 指标 采用 等 值 后 的 阅读 能 力 和 数学 能 力 IRT 分 数 ， 可 以 进 
行 跨 年 级 比较 。 测 验 设计 进一步 的 信息 可 以 参考 相应 的 技术 报告 
数据 中 的 阅读 (9 与 数学 (能 力 的 不 同年 级 的 均值 和 标准 差 以 及 相关 和 矩阵 如 下 所 示 。 


° 


表 2 ECLS-K 学 生 阅读 与 数学 能 力 发 展 均值 、 标 准 差 与 相关 


rl r2 r4 15 тб r7 ml m2 m4 m5 m6 m7 
rl 1.000 
r2 .805 1.000 
r4 .692 .791 1.000 
r5 ‚627 ‚683 .775 1.000 
тб ‚604 654 .732 .862 1.000 
r7 .551 .569 .625 .761 .802 1.000 


ml .784 .735 .689 .673 .657 .597 1.000 

m2 .703 .765 .708 .690 .675 .608 .834 1.000 

m4 598 .658 .729 .691 .677 .606 .732 .796 1.000 

m5 .583 .624 .668 .745 .730 .676 .713 .764 .803 1.000 

m6 .556 .594 .637 .708 .740 .692 .079 28 Лт .881 1.000 

m7 .540 .569 .608 .680 .708 .730 .639 .680 .708 .815 .853 1.000 
M -1.333 -0.754 0.093 0772 1.027 1287 -1.181 -0.696 0041 0.701 1.091 1419 
SD 0526 0512 0463 0.318 0.302 0.392 0.483 0465 0423 0392 0415 0.452 
ТЕ. 相关 系数 均 达 到 p<.001 水 平 。 
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3.2 增长 形态 的 探讨 


首先 ， 探 讨 增长 形态 。 根 据 前 人 研究 (对 源 ， 刘 红 云 ， 


2018)， 该 数据 的 增长 形态 或 存在 关 


键 转折 点 。 故 建立 三 个 模型 进行 比较 : 线性 增长 模型 、 二 次 增长 模型 和 多 阶段 增长 模型 。 将 


能 力 增长 曲线 和 原始 数据 进行 对 比 ， 如 图 10 所 示 。 从 结果 比较 中 可 以 看 出 , 多 阶段 模型 对 于 
数据 的 拟 合 效 果 最 好 (附录 1)， 故 考虑 多 阶段 模型 作为 潜 增 长 形态 。 


2.5 
2 
T 一 一 一 阅读 线性 
1 ORAT > 2 数学 线性 
is Pini 
a o oss sa 2. 
Š ` 一 一 阅读 多 阶段 
一 一 一 一 数学 多 阶段 
х ”阅读 实际 
T o 数学 实际 
-2 


图 10 ECLS-K 学 生 阅 读 、 数 学 能 力 增长 趋势 


3.3 往复 式 影响 的 模型 比较 


在 确定 了 增长 形态 的 基础 上 ,添加 交叉 滞后 部 分 。 建 立 如 表 2 所 示 的 七 个 模型 。 RILCLM 


包含 一 个 随机 截 距 , 可 参照 LGM 建 模 的 方式 , 但 是 注意 
变量 截 距 自 由 估计 (在 LGM 模型 中 一 般 限 定 观测 变量 截 


RI-CLM 限定 潜 变 量 均值 为 0， 观 测 
BA 0, 潜 变量 均值 自由 估计 )。 在 TSE 


模型 中 ， 同 样 不 估计 截 距 因子 的 均值 ， 放 宽 指 标 截 距 自 


估计; 同时， 定义 测量 误差 ， 每 一 


个 指标 都 定义 一 个 潜 变 量 ， 固 定 载荷 为 1; 在 此 基础 上 ， 


自 回 归 和 交叉 滞后 模型 就 可 以 在 淤 


变量 上 进行 定义 了 。ALT 考虑 多 阶段 增长 形态 (转折 点 在 第 三 次 测量 )、 自 回归 与 交叉 滞后 影 


响 : 注意 载荷 矩阵 A 只 影响 第 二 次 到 第 七 次 测量 ， 第 一 次 测量 与 潜 变 量 因子 设 定 成 相关 并 自 
由 估计 。LV-ALT 可 将 TSE 的 语句 与 ALT 的 语句 结合 即 可 (附录 2)。 洲 增长 模型 方面 , LCM-SR 
中 ， 定 义 测 量 误差 为 潜 变 量 并 设 定 “误差 之 误差 ?为 0， 建 立 潜 变 量 自 回 归 和 交叉 清 后 参数 。 


FLCM-SR 中 放宽 “误差 之 误差 "为 0 的 设 定 ( 附 录 3). 
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03 交叉 滞后 潜 增 长 系列 模型 拟 合 (ECLS-K) 
FLCM-SR LCM-SR LV-ALT ` ALT TSE RI-CLM CLM 
AIC 4329299 53602.70 61850.70 60321.31 17715.09 19356.16 35552.48 
BIC 4377027 53976.57 62304.12 60806.54 18112.82 19777.75 35910.44 
-FZ; 26373.52 37616.18 46737.20 45200.76 2258.98 4038.14 18979.92 


ағ 30 43 33 29 40 37 45 

CH 0.841 0.774 0.719 0.704 0.987 0.976 0.975 
RSMEA 0.204 0.204 0.259 0.272 0.051 0.072 0.067 
SRMR 0.203 0.252 1.321 1.309 0.022 0.084 0.059 


ik. FLCM-SR: 因子 结构 化 潜 增 长 模型 ，LCM-SR: 结构 化 潜 增 长 模型 ，LV-ALT: ЖЕН ІШІНІҢ 
长 模型 ，ALT: 自 回归 潜 增 长 模型 ，TSE: 特质 -状态 -误差 模型 ，RI-CLM: 随机 和 截 距 交叉 滞后 模型 ，CLME 
交叉 滞后 模型 . 


n 


七 个 模型 的 模型 拟 合 与 参数 估计 结果 如 表 3 和 表 4 所 示 。 首 先 对 比 CLM、RI-CLM il TSE. 
这 三 个 模型 可 以 看 成 是 逐步 排除 变异 的 过 程 : 添加 稳定 因子 并 排除 测量 误差 。 可 知 添加 了 稳 
定 因子 的 模型 拟 合 优 于 CLM， 而 进一步 排除 测量 误差 之 后 , 模型 拟 合 又 得 到 了 提升 。 然 后 对 
比 ALT 和 LVALT， 后 者 在 前 者 的 基础 上 排除 了 测量 误差 ， 模 型 拟 合 有 所 提升 。LCM-SR 和 
FLCM-SR 的 对 比 和 上 一 对 类 似 ， 后 者 在 排除 了 测量 误差 之 后 ， 模 型 拟 合 得 到 大 幅度 提升 。 
再 次 ， 对 比 添 加 了 增长 趋势 的 四 个 模型 (ELCM-SR、LCM-SR、LVALT 和 ALT)， 整 体 拟 合 都 
不 如 没有 添加 增长 趋势 的 三 个 模型 (TSE、RI-CLM 和 CLM)， 这 也 和 增长 趋势 是 否 完美 拟 合 
数据 有 关 ( 附 录 1)。 最 后 ， 对 比 两 类 基础 模型 的 建 模 (是 否 包含 外 源 变 量 )， 发 现 将 第 一 次 测量 
考虑 成 内 源 变 量 的 模型 (FELCM-SR 和 LCM-SR) 要 比 考虑 成 外 源 变量 的 模型 CVALT 和 ALT) 
拟 合 更 优 。 从 实证 数据 中 综合 来 看 ， 结 合 了 潜 增 长 与 交叉 滞后 的 模型 ， 模 型 设 定 比较 复杂 ， 
参数 限定 较 多 ， 整 体 拟 合 也 较 差 。 整 理 来 看 ， 特 质 -状态 -误差 模型 最 适合 本 数据 的 拟 合 


Жа 交叉 灌 后 潜 增 长 系列 模型 参数 估计 (ECLS-K) 


R2 ON 
КІ 
МІ 


自 回归 交叉 滞后 影响 


FLCM-SR LCM-SR LV-ALT ALT TSE RI-CLM CLM 
Est. S.E. Est. S.E. Est. S.E. Est. S.E. Est. S.E. Est. S.E. Est. S.E. 
1303 0.017 1702 0.042 0.258 0.016 0.353 0042 0.548 0013 042 0.008 0.579 0.007 
-0.391 0.013 -0.919 0.038 0.575 0.015 0.436 0.036 0.349 0.022 0.031 0010 0.286 0.008 
-1.423 0.042 0.531 0.018 -0.660 0.207 0.714 0.013 0.531 0.010 0.403 0.008 0.544 0.007 
1408 0.072 0.462 0.018 1.553 0.213 0.213 0.014 0.296 0.016 -0.003 ™ 0.009 0.250 0.007 
0.347 0.014 -0.670 0.065 -0.978 0.089 0.559 0.017 0.379 0.015 0.169 0.008 0.401 0.006 
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M4 0352 0.010 1157 0.075 1467 0.088 0101 001 0.256 0012 -0.021 0.007 0.203 0.006 


К6 ОМ 
R5 0.889 0010 1.575 0.02 0.318 0.022 0.533 0.038 0.836 0.024 0.252 0.013 0.675 0.007 
M5 0.030 0.009 -0.466 0.014 0.526 0.02 0135 0.023 0.045" 0.016 -0.061 0.009 0.150 0.006 

R7 ON 
R6 0.006" 0.029 1745 0.004 -0.533 0.094 0496 0.071 1004 0.027 0.338 0.017 0.834 0.012 
M6 0658 0.036 -0.400 0.015 1493 0.082 0.394 0.043 0.107 0.022 0.034 0.011 0.206 0.009 

М2 ОМ 
МІ 0274 0.013 1401 0.041 0.633 0.018 0.242 0.077 0.920 0017 0.290 0.010 0.708 0.007 
КІ 1.200 0.049 -0.592 0.043 0.186 0.018 0.572 0.084 -0.042 0.015 0.102 0.007 0.113 0.006 

М4 ОМ 
М2 1046 0.041 0.497 0.017 1.977 0.185 0.840 0.011 0.832 0.025 0.203 0.009 0.647 0.007 
R2 -0.952 0.027 0.342 0013 -1063 0.178 006 0.010 -0.022” 0.019 0.043 0.007 0.096 0.006 


0409 0.011 -0.396 0.037 -0.557 0.095 0.802 0.014 0.806 0.033 0.197 0.008 0.628 0.008 
0.374 0.020 0.661 0.031 1.252 0.096 004” 0.011 0.045" 0.022 0.010 ™ 0.007 0.150 0.007 
-0.049” 0.017 0.905 0.009 0.961 0.020 1.083 0.015 1.105 0.020 0.373 0010 0.840 0.007 
1221 0.023 0.101 0.007 -0.008 ™ 0.024 -0.222 0.014 -0.112 0.021 -0.106 0.012 0.149 0.009 
0.726 0.017 0.884 0.023 1.735 0.074 0.843 0.021 0970 0.013 0.363 0.011 0.792 0.009 
0.029" 0.009 -0.005" 0.003 -0.507 0.089 0.095 0.022 0.073 0.019 -0.046 0.015 0.254 0.012 
a 限制 参数 影响 * 

0.459 0019 0.161 0.026 0371 0024 0578 0007 0.419 001 0322 0.007 0.536 0.003 
0579 0.023 0.349 0014 0.339 0.024 0.074 0.009 0.340 0.018 -0.072 0.006 0.235 0.003 
0210 0.016 0.763 0.034 0.501 0015 0.564 0.013 0.930 0.014 0.214 0.006 0.721 0.003 
0.434 0.016 0.033 0.037 0.196 0.016 0.177 0.009 -0.017” 0.011 -0.044 0.005 0.116 0.003 


ik. ”限制 参数 影响 表示 将 对 应 自 回归 系数 、 交 叉 滞 后 系数 限定 为 相等 ， 估 计 平 均 效 应 ， 即 非 时 变 模 型 。 
ЧЕ DEB АЗ p<.001 水 平 。 


对 比 参 数 估计 ， 不 同 模型 得 出 一 致 性 结论 : 自 回归 系数 比 交 叉 清 后 系数 估计 值 大 。 阅 读 
能 力 的 自 回 归 系 数 比 数学 能 力 的 自 回 归 系 数 略 小 ， 说 明 数 学 能 力 更 倾向 于 惰性 发 展 ， 受 到 先 
前 状态 影响 更 大 。 横 向 对 比 不 同 的 模型 ， 模 型 之 间 的 结果 差异 较 大 。 交 叉 滞 后 参数 部 分 ， 
什 下 结论 。 在 考虑 了 截 距 因 子 消 后 ， 交 叉 滞后 的 参数 佑 
计 值 减 小 ， 而 加 入 了 和 斜率 之 后 ， 自 回归 部 分 略 有 增 大 。 

结合 本 研究 的 实际 情况 ， 在 阅读 和 数学 的 增长 过 程 当中 ， 存 在 着 稳定 的 “能 力 特 质 *， 即 
从 幼儿 园 到 初中 的 发 展 中 保持 个 体 间 的 稳定 性 ， 在 此 基础 之 上 ， 个 体内 发 展 惰性 主要 存在 于 
低 年 级 ， 在 高 年 级 时 个 体内 的 滞留 影响 较 小 ， 容 易 受 到 其 他 能 力 的 影响 ， 阅 读 和 数学 能 力 之 
闻 也 具有 往复 式 影响 。 从 增长 形态 来 说 ， 多 阶段 增长 模型 虽然 比 线性 或 二 次 模型 有 较 好 的 拟 
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合 ， 但 是 不 及 单一 随机 截 距 因子 模型 拟 合 优良 。 综 合 看 来 ， 可 采用 TSE KER. 


4 讨论 与 总 结 


4.1 建 模 思路 : 如 何 分 解 变异 


本 文 探 讨 了 多 元 追踪 研究 中 的 不 同 模型 ， 对 模型 的 构建 和 应 用 进行 了 梳理 ， 并 且 将 模型 
进行 整合 , 拓展 出 具有 操作 化 意义 的 模型 框架 。 结合 以 往 研究 ，Usami 等 人 (2019) 的 框架 并 没 
有 考虑 ALT 模型 的 测量 部 分 , 以 及 LCM-SR 固定 “误差 之 误差 ”为 0 的 假设 。 虽然 Bianconcini 
和 Bolle(2018) 的 模型 考虑 了 ALT 的 测量 误差 ,但 是 目前 的 研究 大 多 以 CLM 为 基础 模型 进行 
ЖЕЛИ, 忽略 了 LGM 为 基础 模型 的 建 模 和 二 者 之 间 的 转换 。 故 本 文 综合 梳理 了 两 个 建 模 思路 ， 
各 自 增 减 相 应 的 变异 成 分 ， 得 到 不 同 的 模型 。 
以 本 文 整合 框架 中 的 FLCM-SR 为 例 ， 首 先 排 除 个 体 间 变异 ， 即 潜 增长 部 分 。 这 部 分 措 
述 了 个 体 随 时 间 发 展 的 稳定 特征 ， 既 描述 了 初始 特质 (基线 )， 又 刻画 了 增长 规律 (成 熟 )。 建 模 
过 程 中 , 潜 增 长 部 分 可 以 被 看 成 是 一 个 需要 被 控制 的 特征 , 由 于 个 体 “ 基 线 ” 和 "成熟 "的 存在 ， 
排除 掉 一 部 分 自然 发 展 的 规律 是 合理 的 (Berry & Willoughby, 2017)。 特别 地 , 在 建立 潜 增 长 模 
| 型 的 时 候 ， 如 果 测 量 的 变量 之 间 没 有 经 过 等 值 ， 则 斜率 因子 的 意义 不 明确 ， 或 只 能 反映 个 体 
> 在 群体 当中 的 相对 变化 (Liu et al., 2020)。 在 本 文 实证 应 用 的 建 模 过 程 中 , BRT CLM 之 外 , 其 
= 实 都 考虑 了 个 体 间 的 成 熟 因 素 ， 只 不 过 RI-CLM 和 TSE 只 加 入 了 截 距 ，ALT/LCM-SR 和 
= LV- ALT/FLCM-SR ЛА. Y Е АЈ ДА ЖЕ ДО ТР НУР R AREE 
[z 的 模型 ， 其 原因 可 能 是 多 阶段 形态 模型 对 本 组 数据 的 拟 合并 非 最 优 的 结果 而 导致。 
在 排除 了 成 熟 影响 之 后 ， 就 可 以 考察 个 体内 的 影响 了 ， 即 FLCM-SR 中 的 “SR( 结 构 化 往 
复式 影响 )”。 在 这 里 ， 个 体内 的 影响 主要 由 自 回 归 和 交叉 溃 后 参数 两 部 分 构成 。 特 别 地 ， 个 
体内 的 自 回归 即 滞留 影响 或 个 体 惰性 ， 研 究 者 应 将 其 和 前 一 段 提 到 的 个 体 间 的 “自然 成 熟 " 进 
行 区 分 (Bery & Willoughby, 2017; Curran & Bauer, 2011; Hamaker et al., 2015; Voelkle et al., 
2014)。“ 谐 留影 响 ” 是 指 个 体内 当前 的 状态 多 大 程度 上 取决 于 他 /她 的 已 有 状态 , 即 变化 的 顽固 
程度 ， 是 个 体 不 变 的 惯性 ， 而 “自然 成 熟 "是 个 体 间 的 发 展 规律 ， 就 像 诸多 心理 学 教材 当中 书 
写 的 “既定 内 容 ” 一 样 。 以 本 文 实证 研究 的 阅读 发 展 为 例 , “自然 成 熟 " 即 潜 增 长 部 分 是 一 般 的 
阅读 发 展 规律 一 一 幼儿 园 到 一 年 级 ， 阅 读 发 展 迅速 ， 进 入 小 学 之 后 ， 阅 读 发 展 速度 降 下 来 。 
而 “滞留 影响 ”是 指 ， 诸 如 一 年 级 的 阅读 能 力 受 到 幼儿 园 阅读 能 力 的 影响 ， 八 年 级 的 阅读 能 力 
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受到 五 年 级 阅读 能 力 的 影响 。 这 种 滞留 影响 在 一 些 灵 活 的 SEM 结构 中 可 以 自由 估计 ， 比 如 
一 年 级 阅读 能 力 受 幼儿 园 阅 读 能 力 的 影响 较 大 ， 而 八 年 级 阅读 能 力 受 五 年 级 阅读 能 力 影 响 较 
小 。 所 以 ， 此 类 模型 区 分 了 个 体内 和 个 体 间 效 应 ， 解 决 了 发 展 当 中 普遍 规律 与 阶段 差异 的 问 
题 ， 更 强调 个 体 在 环境 当中 的 “可 塑性 ”。 

F 面 一 部 分 是 交叉 滞后 影响 ， 即 往复 式 影响 ， 这 部 分 结果 或 许 是 研究 者 最 为 关注 的 焦点 
之 一 。 在 排除 了 个 体 间 变异 之 后 ， 交 叉 清 后 的 影响 或 被 淹没 。 原 因 有 可 能 是 在 比较 传统 的 模 
型 中 ， 交 又 灌 后 的 影响 杂 焰 了 发 展 、 个 体 间 增长 规律 等 因素 ， 使 得 交叉 灌 后 的 影响 被 高 估 
(Berry & Willoughby, 2017)。 但 是 ， 这 种 “高 估 ” 并 非 完全 不 可 接受 ， 因 为 这 种 杂 焰 了 其 他 因素 
的 影响 是 在 某 个 水 平 层面 上 的 解读 ， 可 以 把 它 看 成 是 “整体 影响 ”。 此 类 研究 倾向 于 将 往复 式 
影响 定位 到 个 体内 层面 ， 或 可 以 更 精确 地 表达 如 本 文 一 开始 就 定义 的 “变化 是 如 何 引 起 变化 ” 
的 含义 。 

最 后 的 议题 是 是 否 考虑 测量 误差 ， 即 每 次 观测 变量 的 测验 信和 度 。 本 例 中 LVALT 和 
FLCM-SR 的 拟 合 分 别 略 优 于 ALT 和 LCM-SR， 可 见 考虑 了 测量 信 度 会 一 定 程度 上 提升 模型 
估计 质量 。 值 得 注意 的 是 ， 本 例 中 每 次 测量 都 是 单一 指标 ， 所 以 测量 信 度 是 一 个 固定 值 。 如 
果 每 次 的 测量 有 多 个 指标 ， 或 者 达到 模型 识别 的 要 求 ， 测 量 误差 才 真 正 被 模型 估计 。 可 以 参 
考 Bianconcini 和 Bollen(2018) 提 出 的 LVALT 或 Mulder 和 Hamaker(2020) 提 出 的 多 指标 
КІСІМ 框架 ， 建 构 多 指标 测量 模型 。 综 合 看 来 ， 如 果 有 先 验 信息 ， 加 入 测量 信和 度 的 模型 拟 
合 数据 会 更 精确 ， 但 是 在 只 有 单一 指标 的 情况 下 ， 园 于 模型 识别 等 因素 ， 必 须 人 为 限定 测量 
с 信 度 .此 外 , 在 多 指标 跨 时 测量 时 , 跨 时 测量 等 价 性 也 是 需要 研究 者 提前 考察 的 议题 (Masselink 


= et al., 2018). 


4.2 模型 选择 的 议题 


上 述 论证 的 变异 分 解 理论 建构 应 该 是 研究 者 考虑 的 模型 选择 的 重要 依据 之 一 .尽管 如 此 ， 
模型 选择 的 过 程 中 ， 确 有 一 些 不 稳定 因素 会 影响 研究 者 的 使 用 。 比 如 ， 研 究 者 并 不 知道 " 真 ” 
模型 ， 即 理论 上 真正 的 变异 应 该 如 何 分 解 ， 所 以 无 法 选择 “正确 ”的 模型 。 事 实 上 ， 一 些 模拟 
究 的 结果 也 显示 ， 通 过 不 同 定义 框架 来 生成 的 数据 ， 采 用 “正确 ”模型 的 估计 反 真 率 也 几乎 
只 有 六 七 成 (Usami et al., 2015; 2016)。 这 也 表明 ,在 对 上 述 模型 的 选择 上 ,“ 正 确 ” 模 型 能 被 选 
出 来 的 概率 并 不 是 太 高 ， 而 “错误 ”模型 使 用 之 后 ， 其 实 对 结果 解释 并 不 会 造成 太 大 的 影响 。 
所 以 ， 在 实证 研究 中 ， 并 不 是 模型 越 复杂 拟 合 越 佳 ， 系 列 的 拟 合 指标 是 研究 者 选择 模型 的 首 
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要 参考 依据 (Usamietal., 2015). 

建立 复杂 的 模型 容易 ， 但 是 选择 合适 的 模型 并 不 容易 。 研 究 者 应 首先 根据 理论 基础 设 定 
合理 的 模型 。 如果 某 个 特质 因素 确 是 稳定 特质 , 则 考虑 包含 随机 截 距 的 模型 以 提出 稳定 特质 ; 
如 果 测 量 时 间 较 长 且 确 实 有 发 展 (如 能 力 )， 再 考虑 增长 因子 。 一 般 而 言 ， 提 取出 稳定 特质 的 
模型 会 具有 较 好 的 模型 拟 合 ， 而 增长 因子 与 交叉 滞后 部 分 夹杂 在 一 起 会 影响 带 给 研究 者 一 些 
解释 困难 。 其 次 ， 通 过 备 选 模型 之 间 的 拟 合 指标 比较 来 选择 模型 是 应 用 研究 者 应 该 参考 的 主 
要 评判 标准 之 一 。 此 类 模型 由 于 复杂 程度 差异 较 大 ， 前 人 模拟 研究 结果 建议 更 多 采用 CFI 和 
RMSEA 来 选择 模型 ， 信 息 函数 的 表现 并 不 理想 ， 且 BIC 容易 惩罚 复杂 模型 而 选择 更 简单 的 
模型 (Usami et aL, 2015)。 综合 以 往 研究 的 建议 Berry & Willoughby, 2017; Usami et al., 2019) 和 
本 文 实证 数据 的 结果 ， 应 用 者 应 使 用 诸如 CFI. ТІЛ. RMSEA, SRMR 以 及 AIC 等 指标 进行 
综合 判断 。 


Қа! 


4.2 其 他 议题 与 研究 展望 


对 于 模型 的 选择 ， 近 期 有 学 者 提出 “结构 方程 模型 树 (SEM tree)”" 和 “结构 方程 模型 森林 
(SEM forest)” 借鉴 了 管理 心理 学 当中 的 决策 树 概 念 , 来 对 模型 的 选择 提供 证 据 (Brandmaier et 
o al., 2016; Brandmaier, Von et al., 2013)。 未 来 可 以 通过 模拟 研究 来 给 出 模型 选择 的 建议 ， 比 如 
I 如 何 选择 真 模型 ， 模 型 误 设 的 代价 ， 以 及 个 体 间 /个 体内 差异 达到 何 种 程度 的 时 候 推荐 哪 类 模 
型 等 。 

另 一 方面 ， 时 变 模 型 的 构建 可 以 使 得 模型 估计 更 为 精确 。 比 如 自 回归 模型 中 的 非 时 变 影 
€) 响 在 SEM 框架 下 可 以 将 参数 进行 一 定 程度 的 放宽 ; 滞留 影响 所 在 不 同 测量 次 数 之 间 可 以 自 
由 估计 ， 模 型 限定 更 少 ， 可 以 有 效 提高 模型 整体 拟 合 。Bringmann 等 人 (2018) 提 出 了 时 变 自 回 
归 模 型 (time-variant autoregressive model), ， 假 设 自 回归 影响 在 不 同 次 的 时 间 间 隔 不 同 
(Bringmann et al., 2018)。 若 设 定 滞留 效应 的 非 时 变性 的 理论 假设 并 不 确定 ， 则 可 以 放宽 该 假 
设 ， 整 体 拟 合 更 优 (Dinget al., 2020)。 但 是 需要 注意 参数 识别 的 问题 。 

此 外 ,研究 者 在 追踪 研究 中 根据 不 同 测量 时 间 间 隔 ， 对 于 时 间 变 量 的 精度 要 求 会 有 所 不 
FJ. ТЕ SEM HEAR F, BIERE A 大 多 都 固定 测量 时 间 。 这 会 增 大 估计 误差 ， 降 低 模型 的 敏 
感度 (Liu et al., 2015; Sterba, 2014)。 在 此 基础 上 , 借鉴 多 水 平 模型 框架 ， 可 以 放宽 测量 时 间 固 


定 的 假设 ，SEM 也 逐渐 开始 使 用 连续 时 间 模 型 (continuous time models，CTM)， 将 被 试 参 与 


的 时 间 定 义 为 一 个 随机 数 ， 而 非 固 定 的 测量 次 数 (Deboeck & Preacher, 2016; Hamaker & 
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Wichers, 2017; Schuurman et al., 2016). 在 新 版 本 的 SEM 软件 (如 Mplus 8.0), і 


经 得 到 技术 支撑 ， 可 以 供 研究 者 灵活 选用 。 
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Abstract: When conducting the multivariate longitudinal studies, reciprocal relationship and 
latent trajectory are two of the focusing issues. These two issues could be flexibly combined by 
other research questions, such as the measurement error, the random factor, as well as the 
combination of the above topics. Such a combination yields a more complex model definition 
exploring the longitudinal relations, such as factor cross-lagged model, random-intercept 
cross-lagged model, trait-state-error model, autoregressive trajectory model, etc. In the study, a 
factor latent curve model with structured reciprocals model was built as an extension and unified 
framework including all the components discussed above. The empirical dataset, Early Childhood 
Longitudinal Survey-Kindergarten (ECLS-K), was used as an illustrating example. Results 
indicated that the trait-state-error model best described the data. Finally, we summarized how the 
results could be interpreted and offered suggestions on model selection for the researchers. 

Key words: longitudinal study, reciprocal effect, growth trajectory, factor latent curve model 
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附录 1 线性 、 二 次 、 多 阶段 增长 模型 结果 参数 估计 


线性 增长 模型 (LGM) 


Est. 


Means 


INTR -0.889 


SLPR 0.353 


QUAR 


SLPR2 


INT™ -0.865 


SLPM 0.363 


QUAM 


SLPM2 


Variances 


INTR 0.205 


SLPR 0.004 


QUAR 


INT™ 0.184 


SLPM 0.003 
QUAM 

SLPM2 

Residual 

rl 0.288 
I2 0.065 
r4 0.266 
I9 0.202 
тб -0.009 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


S.E. 


0.004 


0.001 


0.003 


0.001 


0.003 


0.002 


0.004 


0.002 


0.003 


0.003 


0.001 
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二 次 增长 模型 (QGMD 
Est. S.E. 
-L112 *** 0.004 
0.720 *** 0,002 
-0,055 с 0 
-1.024 *** 0.003 
0.64 *** 0.001 
-0.042 с 0 
0.244 *** 0,003 
0.015 *** 0.001 
0 жж 0 
0.196 *** 0.002 
0.008  *** 0 
0 жж 0 
0.131 *** 0,002 
0.041 *** 0.001 
012 *** 0.002 
0.014 *** 0001 
0.060  *** 0.001 


Est. 


-1.321 


1.158 


-1.009 


-1.167 


0.914 


-0.703 


0.24 


0.026 


0.014 


0.2 


0.008 


-0.003 


0.045 


0.059 


0.18 


0.026 


0.008 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


жж 


多 阶段 增长 模型 PGMD) 


S.E. 


0.004 


0.003 


0.004 


0.003 


0.002 


0.003 


0.003 


0.001 


0.002 


0.002 


0.001 


0.001 


r7 0.657 *** 0.01 0.18 *** 0.008 0.00 *** 0.002 


ml 0.163 *** 0002 0.083  *** 0.001 0.042  *** 0.001 
m2 0.035  *** 0.001 0.029 *** (0001 0.0385  *** 0.001 
m4 0.179 *** 0.002 0.088 *** 0.001 0.082  *** 0.001 
m5 012 *** 0002 0017 *** (0001 003 *** 0.001 
m6 0.002 0.001 0.041 *** 0.001 0.024 *** 0.001 
m7 0.617 *** 0.01 0.045  *** 0.004 0.134  *** 0.003 
SLPR WITH 

INTR -0.02  *** 0 -0.052 *** 0.001 -0.069 *** 0.002 
QUAR 0.004 *** 0 

SLPR2 0.06 *** 0.002 
SLPR WITH 

QUAR -0.001  *** 0 

SLPR2 -0.02 *** 0.001 
SLPM WITH 

SLPM -0.01  *** 0 -0.023  *** 0.001 -0.031 *** 0.001 
ОСАМ 0.002. *** 0 

SLPM2 0.026 *** 0.001 
SLPM WITH 

QUAM -0.001 *** 0 

SLPM2 -0.003 *** 0.001 

жж p<.001 


附录 2 潜 变量 自 回归 潜 增长 模型 Mplus 语句 与 示意 图 


TITLE: this is an example of LV-ALT 
DATA: FILE IS eclsk theta-alt.dat; 
VARIABLE: NAMES ARE 
CIRARTHT R CIRAMTHT R CIRGTHT R C2RARTHT R C2RAMTHT R C2RGTHT R C3RARTHT R 


C3RAMTHT К C3RGTHT К CARARTHT К CARAMTHT К САКСТНТ К CSRARTHT К CSRAMTHT К 
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CSR2STHT К C6R4RTHT_R СӨКАМТНТ R C6R2STHT К C7R4RTHT_R C7RAMTHT R C7R2STHT R 


PILEARN P2LEARN P4LEARN TILEARN T2LEARN TALEARN TSLEARN TóLEARN; 


USEVARIABLES ARE 


CIR4RTHT_R C2R4RTHT R CARARTHT К C5R4RTHT_R C6R4RTHT_R C7RARTHT К 


CIRAMTHT К C2RAMTHT К C4R4MTHT_R C5R4MTHT_R C6R4MTHT_R C7R4MTHT R; 


MISSING AREALL (-9); 


MISSING AREALL (-8); 


MISSING AREALL (-7); 


MISSING AREALL (-1); 

ANALYSIS: MODEL = NOCOVARIANCES; 
estimator=mlr; 

miteration = 5000;convergence=0.01; 


MODEL: 


Idefine the growth( 排 除 第 一 次 测量 为 外 源 变 量 后 ， 定 义 多 阶段 潜 增长 模型 ， 前 三 次 测量 为 斜率 1, 26 


后 为 斜率 2); 
INTR by C2R4RTHT_R@1 C4R4RTHT_R@1 С5КАКТНТ К@ 1 CORARTHT К@ І CTRARTHT RG1; 
SLPR by C2R4RTHT_R@0.5 C4R4RTHT_R@1.5 
CSRARTHT К@3.5 C6R4RTHT_R@5.5 CTRARTHT R 8.5; 
SLPR2 Бу C2RARTHT К@0 CARARTHT RG0 
CSRARTHT К@2 CGRARTHT R24 СТКАКТНТ К @7; 
INTM by C2R4MTHT_R@1 CARAMTHT R1 CSRAMTHT RG 1 C6R4MTHT_R@1 СТЕАМТНТ КІ; 
SLPM by C2R4MTHT_R@0.5 САКАМТНТ К@ 1.5 
CSRAMTHT RG3.5 СӨКАМТНТ К@5.5 CTRAMTHT RG8.5; 
SLPM2 by C2RAMTHT К @0 CARAMTHT RG0 
С5КАМТНТ RG2 C6RAMTHT R4 СТЕАМТНТ RG7; 
[C2RARTHT R-C7RAMTHT R0]; 
INTR SLPR INTM SLPM slpr2 slpm2; 
[INTR SLPR INTM SLPM slpr2 slpm2]; 
intr with slpr slpr2; slpr with slpr2; 


intm with slpm slpm2; slpm with slpm2; 
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!define state factor( 定 义 状 态 变量 ， 即 每 次 测量 的 信 度 ); 


R1 BY CIR4RTHT_R@1; R2 BY CORARTHT RG1; КА BY CARARTHT К@ 1; 


R5 BY CSRARTHT RG1; R6 BY CGORARTHT RG1; R7 BY CTRARTHT R1; 


МІ BY CIRAMTHT RG1; M2 BY C2RAMTHT RG 1; МА BY САКАМТНТ КӨЗІ; 


M5 BY С5ЕАМТНТ RG1; M6 BY СӨКАМТНТ RG 1; M7 BY СТЕАМТНТ КІ; 


>l 


[El 


!define the exogenous indicators( 定 义 第 一 次 测量 的 因子 为 外 源 变量 ); 
R1 with intr slpr slpr2; 
MI with intm slpm slpm2; 


E 


!define measurement error( 定 义 误差 变量 ); 


CIRARTHT R-C7RARTHT R(l); 


CIRAMTHT R-C7RAMTHT R(2); 


!define auto-regression and cross-lagged regression( 定 义 自 回归 与 交叉 滞后 影响 ); 
R20NRI; R4 ON R2; R5 ON R4; 
R6 ON R5; R7 ON R6; 


M2 ON MI; M4 ON M2; М5 ON M4; 


M6 ON M5; M7 ON M6; 


M2 ON RI; M4 ON R2; M5 ON R4; 


M6 ON R5; M7 ON R6; 
R2 ON M1; R4 ON M2; R5 ОМ M4; 
R6 ON M5; R7 ON M6; 


!define the residuals (4E COWRA 77 28. PINAR); 


R1; R2-R7(3); 


М1; M2-M7(4); 


!define the correlation between residuals( 定 义 因子 之 间 协 方差 /相关 ); 


КІ WITH МІ; R2 WITH М2; R4 WITH М4; 


R5 WITH M5; R6 WITH M6; R7 WITH M7; 


OUTPUT: STDYX;tech4;tech1; 
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附录 3 因子 结构 化 潜 增 长 模型 Mplus 语句 与 示意 图 


TITLE: this is an example of FLCM-SR 

DATA: FILE IS eclsk_theta+alt.dat; 

VARIABLE: NAMES ARE 
CIRARTHT R CIRAMTHT R CIRGTHT R C2RARTHT R C2RAMTHT R C2RGTHT R C3RARTHT R 
C3RAMTHT R C3RGTHT R САКАКТНТ К CARAMTHT К C4RGTHT_R C5R4RTHT R С5КАМТНТ R 
CSR2STHT R C6R4RTHT_R C6RAMTHT R C6R2STHT R C7RARTHT К C7RAMTHT R C7R2STHT R 
PILEARN P2LEARN P4LEARN TILEARN T2LEARN T4LEARN TSLEARN T6LEARN; 

USEVARIABLES ARE 
CIRARTHT К C2RARTHT R CARARTHT К CSRARTHT К CGRARTHT R C7RARTHT К 
CIRAMTHT К C2RAMTHT К C4RAMTHT К CSRAMTHT К C6RAMTHT Е C7RAMTHT R; 

MISSING AREALL (-9); 

MISSING AREALL (-8); 


MISSING AREALL (-7); 


MISSING AREALL (-1); 
ANALYSIS: MODEL = NOCOVARIANCES; 
estimator-mlr; 


miteration = 5000;convergence=0.01; 
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MODEL: 


!define the growth( 定 义 多 阶段 潜 增 长 模型 ， 前 三 次 测量 为 斜率 1， 第 三 次 以 后 为 斜率 2); 


ІШ 


INTR by CIRARTHT К@І C2R4RTHT_R@1 C4R4RTHT_R@1 CSR4RTHT_R@1 C6R4RTHT_R@1 
C7RARTHT К@ 1; 
SLPR by CIRARTHT К@0 C2RARTHT К @0.5 CARARTHT К@ 1.5 
С5КАКТНТ К@3.5 СӨКАКТНТ К@ 5.5 СТЕАКТНТ К(8.5; 
SLPR2 by CIRARTHT К@0 C2R4RTHT_R@0 C4R4RTHT_R@0 
C5R4RTHT_R@2 C6R4RTHT_R@4 C7R4RTHT_R@7; 
INTM by CIR4MTHT_R@1 C2R4MTHT_R@1 C4R4MTHT_R@1 CSR4MTHT_R@1 C6R4MTHT_R@1 
C7R4MTHT_R@1; 
SLPM by CIR4MTHT_R@0 C2R4MTHT_R@0.5 C4R4MTHT_R@1.5 
CSRAMTHT RG3.5 CORAMTHT К@5.5 C7R4MTHT_R@8.5; 
SLPM2 by CIRAMTHT R0 C2RAMTHT R0 CARAMTHT RG0 
С5КАМТНТ RG2 C6R4MTHT_R@4 СТЕАМТНТ RG7; 
[CIRARTHT R-C7RAMTHT R0]; 
INTR SLPR INTM SLPM slpr2 slpm2; 
[INTR SLPR INTM SLPM slpr2 slpm2]; 
intr with slpr slpr2; slpr with slpr2; 


intm with slpm slpm2; slpm with slpm2; 


!define state factor( 定 义 状 态 变量 ， 即 每 次 测量 的 信 度 ); 


R1 BY CIRARTHT RG1; R2 BY CORARTHT RG1; КА BY CARARTHT К@ 1; 


R5 BY CSRARTHT RG1; R6 BY CGORARTHT RG1; R7 BY CTRARTHT R1; 


МІ BY CIRAMTHT RG1; M2 BY C2RAMTHT R91; МА BY CARAMTHT К@ 1; 


M5 BY C5R4AMTHT_R@1; M6 BY СӨКАМТНТ RG 1; M7 BY СТЕАМТНТ КІ; 


!define measurement error( 定 义 误差 变量 ); 
CIRARTHT R-C7RARTHT R; 


CIRAMTHT R-C7RAMTHT R; 


ІШ 


要 与 交叉 滞后 影响 ); 


!define auto-regression and cross-lagged regression(XE X. Fl [Hl 
R2 ON R1; R4 ON R2; R5 ON R4; 


R6 ON R5; R7 ON R6; 
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M2 ON MI; M4 ON M2; М5 ON М4; 
M6 ON M5; M7 ON M6; 

M2 ON RI; M4 ON R2; М5 ON R4; 
M6 ON R5; M7 ON R6; 

R2 ON M1; RA ON M2; R5 ON M4; 

R6 ON M5; R7 ON M6; 

!define the residuals (定义 被 试 内 残 差 ); 
R1; К2-К7(3); 


М1; M2-M7(4); 


!define the correlation between residuals( 定 义 因 子 之 间 协 方差 /相关 ); 
R1 WITH МІ; R2 WITH М2; R4 WITH M4; 
R5 WITH M5; R6 WITH M6; R7 WITH M7; 


OUTPUT: STDYX;tech4;tech1; 
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